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Abstract.  The  problem  of  tropical  cyclone  formation  re¬ 
quires  among  other  things  an  improved  understanding  of 
recirculating  flow  regions  on  sub- synoptic  scales  in  a  time 
evolving  flow  with  typically  sparse  real-time  data.  This  recir¬ 
culation  problem  has  previously  been  approached  assuming 
as  a  first  approximation  both  a  layer-wise  two-dimensional 
and  nearly  steady  flow  in  a  co-moving  frame  with  the 
parent  tropical  wave  or  disturbance.  This  paper  provides 
an  introduction  of  Lagrangian  techniques  for  locating  flow 
boundaries  that  encompass  regions  of  recirculation  in  time- 
dependent  flows  that  relax  the  steady  flow  approximation. 

Lagrangian  methods  detect  recirculating  regions  from 
time-dependent  data  and  offer  a  more  complete  methodology 
than  the  approximate  steady  framework.  The  Lagrangian  ref¬ 
erence  frame  follows  particle  trajectories  so  that  flow  bound¬ 
aries  which  constrain  particle  transport  can  be  viewed  in 
a  frame-independent  setting.  Finite-time  Lagrangian  scalar 
field  methods  from  dynamical  systems  theory  offer  a  way 
to  compute  boundaries  from  grids  of  particles  seeded  in  and 
near  a  disturbance. 

The  methods  are  applied  to  both  a  developing  and  non¬ 
developing  disturbance  observed  during  the  recent  pre¬ 
depression  investigation  of  cloud  systems  in  the  tropics 
(PREDICT)  experiment.  The  data  for  this  analysis  is  de¬ 
rived  from  global  forecast  model  output  that  assimilated  the 
dropsonde  observations  as  they  were  being  collected  by  re¬ 
search  aircraft.  Since  Lagrangian  methods  require  trajectory 
integrations,  we  address  some  practical  issues  of  using  La¬ 
grangian  methods  in  the  tropical  cyclogenesis  problem.  La¬ 
grangian  diagnostics  are  used  to  evaluate  the  previously  hy¬ 
pothesized  import  of  dry  air  into  ex-Gaston,  which  did  not 
re-develop  into  a  tropical  cyclone,  and  the  exclusion  of  dry 
air  from  pre-Karl,  which  did  become  a  tropical  cyclone  and 
later  a  major  hurricane. 


1  Introduction 

1.1  The  genesis  problem  and  pouch  hypothesis 

The  prediction  of  tropical  cyclogenesis1  is  an  important 
problem  in  tropical  meteorology  and  geophysical  fluid  dy¬ 
namics.  While  much  effort  has  been  devoted  to  the  study  of 
mature  tropical  cyclones,  and  their  track  and  intensity  esti¬ 
mates,  only  recently  has  a  significant  research  effort  been 
dedicated  to  the  problem  of  the  transformation  of  a  tropical 
disturbance  into  a  tropical  depression. 

A  new  paradigm  for  tropical  cyclogenesis,  dubbed  the 
“marsupial  paradigm”,  was  recently  developed  by  Dunker- 
ton,  Montgomery,  and  Wang  (Dunkerton  et  al.,  2009)  (here¬ 
after  DMW09).  This  work,  and  that  of  Wang  (2012)  showed 
that  a  region  of  recirculation  along  a  propagating  Easterly 
wave  occurred  in  all  pre-genesis  cases  analyzed  therein.  The 
center  of  the  sub- synoptic  scale  circulation,  defined  as  the 
intersection  of  the  parent  wave’s  trough  axis  and  the  wave’s 
critical  latitude  in  the  lower  troposphere,  was  demonstrated 
to  be  the  preferred  location  of  cyclogenesis.  The  region  of 
recirculation,  coined  the  “pouch”,  protects  the  low-  to  mid¬ 
level  embryonic  vortex  from  adverse  conditions,  while  the 

^he  glossary  on  NOAA’s  Hurricane  Research  Division’s  web¬ 
site  uses  tropical  cyclone  as  the  generic  term  for  a  nonfrontal 
synoptic- scale  low-pressure  system  over  tropical  or  sub- tropical  wa¬ 
ters  with  organized  convection  (i.e.  thunderstorm  activity)  and  a 
definite  cyclonic  surface  wind  circulation.”  Notably,  this  definition 
does  not  invoke  any  wind  threshold.  The  same  glossary  defines  a 
“tropical  depression”  as  a  tropical  cyclone  with  maximum  sustained 
surface  winds  of  less  than  17  ms-1  (34  kt,  39  mph)  and,  in  the  At¬ 
lantic  and  Eastern  Pacific  Basins,  a  “tropical  storm”  as  a  tropical 
cyclone  with  surface  winds  between  17  m  s_  1  and  33  m  s“ 1 .  In  this 
study  we  will  define  genesis  as  the  formation  of  a  tropical  depres¬ 
sion  and  we  impose  no  formal  threshold  on  wind  speed. 


Published  by  Copernicus  Publications  on  behalf  of  the  European  Geosciences  Union. 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

03  DEC  2012 


2.  REPORT  TYPE 


3.  DATES  COVERED 


00-00-2012  to  00-00-2012 


4.  TITLE  AND  SUBTITLE 


A  Lagrangian  analysis  of  a  developing  and  non-developing  disturbance 
observed  during  the  PREDICT  experiment 


6.  AUTHOR(S) 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Postgraduate  School  (NPS), Monterey, CA, 93943 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


5f.  WORK  UNIT  NUMBER 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

The  problem  of  tropical  cyclone  formation  requires  among  other  things  an  improved  understanding  of 
recirculating  flow  regions  on  sub-synoptic  scales  in  a  time  evolving  flow  with  typically  sparse  real-time 
data.  This  recirculation  problem  has  previously  been  approached  assuming  as  a  first  approximation  both  a 
layer-wise  two-dimensional  and  nearly  steady  flow  in  a  co-moving  frame  with  the  parent  tropical  wave  or 
disturbance.  This  paper  provides  an  introduction  of  Lagrangian  techniques  for  locating  flow  boundaries 
that  encompass  regions  of  recirculation  in  timedependent  flows  that  relax  the  steady  flow  approximation. 
Lagrangian  methods  detect  recirculating  regions  from  time-dependent  data  and  offer  a  more  complete 
methodology  than  the  approximate  steady  framework.  The  Lagrangian  reference  frame  follows  particle 
trajectories  so  that  flow  boundaries  which  constrain  particle  transport  can  be  viewed  in  a 
frame-independent  setting.  Finite-time  Lagrangian  scalar  field  methods  from  dynamical  systems  theory 
offer  a  way  to  compute  boundaries  from  grids  of  particles  seeded  in  and  near  a  disturbance.  The  methods 
are  applied  to  both  a  developing  and  nondeveloping  disturbance  observed  during  the  recent  predepression 
investigation  of  cloud  systems  in  the  tropics  (PREDICT)  experiment.  The  data  for  this  analysis  is  derived 
from  global  forecast  model  output  that  assimilated  the  dropsonde  observations  as  they  were  being  collected 
by  research  aircraft.  Since  Lagrangian  methods  require  trajectory  integrations,  we  address  some  practical 
issues  of  using  Lagrangian  methods  in  the  tropical  cyclogenesis  problem.  Lagrangian  diagnostics  are  used 
to  evaluate  the  previously  hypothesized  import  of  dry  air  into  ex-Gaston,  which  did  not  re-develop  into  a 
tropical  cyclone,  and  the  exclusion  of  dry  air  from  pre-Karl,  which  did  become  a  tropical  cyclone  and  later 
a  major  hurricane. 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 


17.  LIMITATION  OF 
ABSTRACT 


19a.  NAME  OF 
RESPONSIBLE  PERSON 


18.  NUMBER 
OF  PAGES 


unclassified 


b.  ABSTRACT 

unclassified 


c.  THIS  PAGE 

unclassified 


Same  as 
Report  (SAR) 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


11356 


B.  Rutherford  and  M.  T.  Montgomery:  A  Lagrangian  analysis  during  PREDICT 


proto- vortex  within  is  nurtured  and  amplified  by  convection 
within  the  pouch.  The  three  main  hypotheses  of  the  marsu¬ 
pial  paradigm  can  be  summarized  as  follows: 

1.  The  cat’s  eye  region  within  the  critical  layer  has  cy¬ 
clonic  recirculation  with  little  strain  deformation,  and 
its  center  serves  as  the  focal  point  for  vorticity  aggrega¬ 
tion  of  convectively  generated  cyclonic  vertical  vorticity 
and  their  vortical  remnants. 

2.  Air  inside  the  cat’s  eye  is  repeatedly  moistened  by  con¬ 
vection  and  is  protected  to  some  extent  from  the  lateral 
intrusion  of  dry  air. 

3.  The  parent  wave  is  maintained  and  possibly  enhanced 
by  convectively  amplified  mesoscale  eddies  within  the 
wave. 

1.2  Cases  from  the  PREDICT  experiment 

The  marsupial  paradigm  was  tested  in  part  during  the  Pre- 
Depression  Investigation  of  Cloud  systems  in  the  Tropics 
(PREDICT)  experiment.  The  PREDICT  experiment  was  de¬ 
signed  to  investigate  the  problem  of  genesis  during  the  2010 
Atlantic  hurricane  season  through  a  series  of  flights  into  and 
data  analysis  of  developing  and  non-developing  tropical  dis¬ 
turbances.  Along  with  flight-level  and  dropsonde  data,  ob¬ 
servational  satellite  data  was  analyzed  and  an  ensemble  of 
numerical  weather  forecast  models  were  run  each  day  during 
the  experiment.  A  description  of  the  overall  experiment,  as 
well  as  a  summary  of  some  first  results,  is  given  by  Mont¬ 
gomery  et  al.  (2012). 

In  a  follow-up  paper,  Smith  and  Montgomery  (2012)  car¬ 
ried  out  a  more  in-depth  analysis  of  the  convective  environ¬ 
ment  in  a  subset  of  developing  and  non-developing  tropi¬ 
cal  disturbances  during  the  last  four  week  period  of  the  ex¬ 
periment.  In  all  systems,  there  is  a  tendency  for  the  lower 
troposphere  to  moisten,  but  in  the  system  that  did  not  de¬ 
velop,  the  middle  and  upper  troposphere  became  progres¬ 
sively  drier  during  the  five  missions.  The  most  prominent  dif¬ 
ference  between  the  non-developing  system  and  the  two  sys¬ 
tems  that  developed  was  the  much  larger  reduction  of  equiva¬ 
lent  potential  temperature  between  the  surface  and  a  height  of 
3  km,  typically  25  K  in  the  non-developing  system,  compared 
with  only  17  K  in  the  systems  that  developed.  Conventional 
wisdom  would  suggest  that,  for  this  reason,  the  convective 
downdraughts  would  be  stronger  in  the  non-developing  sys¬ 
tem  and  would  thereby  act  to  suppress  the  development. 
Smith  and  Montgomery  (op.  cit.)  presented  an  alternative  hy¬ 
pothesis  that  the  drier  mid-level  air  weakens  the  convective 
updraughts  and  thereby  weakens  the  amplification  of  sys¬ 
tem  relative  vorticity  necessary  for  development.  Although 
detailed  testing  of  the  alternative  hypothesis  is  forthcoming 
(Kilroy  and  Smith,  2012),  the  findings  of  Montgomery  et 
al.  (2012)  and  Smith  and  Montgomery  (2012)  raise  the  im¬ 
portant  question  of  what  caused  the  progressive  drying  for 


the  non-developing  disturbance  and  what  prevented  drying 
in  the  developing  disturbances?  As  argued  in  DMW09,  the 
answer  to  this  question  requires  the  realization  of  the  La¬ 
grangian  nature  of  the  genesis  problem  and  the  importance 
of  identifying  intrinsic  flow  boundaries  that  limit/enable  the 
transport  and  mixing  of  adiabatic  invariants  between  the  po¬ 
tential  development  region  and  its  environment. 

The  role  of  structures  and  boundaries  is  important  for  un¬ 
derstanding  the  transport  of  moisture  and  absolute  vertical 
vorticity,  two  important  elements  in  the  genesis  of  a  trop¬ 
ical  cyclone  (DMW09).  Considering  this  problem  in  time- 
dependent  velocity  fields  could  ultimately  lead  to  better  pre¬ 
diction  of  cyclogenesis. 

The  cases  that  are  analysed  in  this  study  are  the  developing 
case  of  Karl  and  the  nondeveloping  case  of  Gaston.  The  over¬ 
all  evolution  of  the  disturbances  can  be  seen  in  Fig.  4,  which 
shows  the  evolution  of  relative  humidity,  the  Okubo- Weiss 
parameter,  and  vorticity  averaged  over  a  3  degree  radius  cir¬ 
cle  surrounding  the  circulation  center.  Karl  remained  moist  at 
all  levels  while  Gaston  started  to  dry  at  700  hPa  and  above. 
One  of  the  outcomes  of  this  work  will  be  an  explanation  of 
the  drying  in  Gaston  and  the  moistening  in  the  Karl  distur¬ 
bance  due  to  the  interaction  of  the  disturbances  with  their 
environments  through  the  inspection  of  flow  boundaries. 

The  study  of  flow  kinematics  during  genesis  is  important 
for  validating  both  aspects  of  the  second  hypothesis  of  the 
marsupial  paradigm  for  candidate  disturbances.  The  second 
hypothesis  requires  Lagrangian  recirculation  along  the  trav¬ 
elling  parent  wave  to  protect  the  pouch  and  also  requires  that 
the  lateral  intrusion  of  dry  air  be  minimal.  Both  of  these  cri¬ 
teria  can  be  verified  by  finding  persistent  boundaries  at  the 
edge  of  the  primary  circulation  which  allow  recirculation 
within  the  nascent  vortex  and  limit  interaction  from  the  envi¬ 
ronment.  Recirculation  of  fluid  parcels  in  a  time-dependent 
flow  does  not  necessarily  occur  in  regions  of  instantaneous 
positive  cyclonic  vorticity,  but  only  occurs  when  the  parcel 
remains  entrained  for  a  sufficiently  long  time. 

The  Okubo-Weiss  (OW)  criterion  (as  defined  herein)  and 
vorticity  are  Eulerian  measures  of  rotation  commonly  used 
in  atmospheric  predictions  and  analysis.  Vorticity  is  defined 
mathematically  as  the  curl  of  the  velocity  field  and  represents 
twice  the  local  rotation  rate  of  a  fluid  particle.  The  (OW)  cri¬ 
terion  helps  characterize  vortex  resilience  by  measuring  the 
relative  role  of  vorticity  in  relation  to  strain  and  shear,  and 
helps  reveal  particles  that  have  high  vorticity  and  low  strain 
deformation  and  remain  entrained  in  vortices. 

1.3  Flow  boundaries 

Flow  boundaries  in  the  Eulerian  or  Earth-relative  frame  of 
reference  are  seen  by  streamlines.  These  streamlines  are 
commonly  used  by  forecasters  to  declare  a  tropical  depres¬ 
sion  when  closed  circulation  is  seen.  However,  it  is  well 
known  that  particle  trajectories  may  cross  Eulerian  stream¬ 
lines  in  time  dependent  flows.  Developing  disturbances  can 
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alternatively  be  viewed  in  a  storm-relative  frame  of  refer¬ 
ence  travelling  at  the  propagation  speed  of  the  parent  East¬ 
erly  wave.  We  refer  to  this  reference  frame  as  the  co-moving 
frame.  In  this  co-moving  frame,  streamlines  often  show 
closed  circulation  in  the  wave  trough  at  sufficiently  large 
wave  amplitude  prior  to  the  declaration  of  a  tropical  depres¬ 
sion.  The  use  of  the  co-moving  frame  has  led  to  a  clearer 
picture  of  the  background  flow  geometry  that  favors  subse¬ 
quent  development  of  a  disturbance.  However,  this  viewpoint 
is  still  not  widely  used  in  forecasting  applications. 

Streamlines  in  a  co-moving  reference  frame  formed  flow 
boundaries  for  a  steady  tropical  cyclone  vortex  studied  by 
Riemer  and  Montgomery  (2011),  and  were  shown  to  influ¬ 
ence  the  import  of  dry  air  to  the  inner-core  of  the  vortex. 
Similar  boundaries  were  found  during  the  PREDICT  experi¬ 
ment  by  Evans  et  al.  (201 1). 

In  the  cases  studied  during  PREDICT,  the  time- variation 
of  the  translation  speed  was  sometimes  sufficiently  large  that 
streamlines  did  not  exhibit  consistently  closed  circulation  un¬ 
der  small  variations  in  the  translation  speed.  In  addition,  the 
co-moving  frame  does  not  account  for  time-dependence  of 
velocities,  so  the  pouch  boundary  in  the  co-moving  frame 
may  not  be  impermeable  to  the  transport  of  air  parcels.  The 
deficiencies  of  the  co-moving  frame  can  be  resolved  by  con¬ 
sidering  the  Lagrangian  frame  of  reference  following  particle 
trajectories. 

In  contrast  to  the  co-moving  frame  following  a  feature 
which  assumes  a  feature-relative  steady  flow,  the  Lagrangian 
reference  frame  can  be  used  to  locate  persistent  flow  bound¬ 
aries  in  flows  with  time-dependence. 

If  the  flow  is  steady,  the  stable  and  unstable  manifolds  of 
a  hyperbolic  fixed  point,  or  saddle-point  partition  the  flow. 
The  manifolds  are  distinguished  from  other  streamlines  since 
they  show  maximal  linear  instability  for  particles  located 
close  to  the  manifolds. 

The  time-dependent  analog  of  a  hyperbolic  fixed  point  is 
called  a  hyperbolic  trajectory  and  is  often  associated  with 
persistent  though  not  necessarily  fixed  stagnation  points  of 
the  steady  flow.  The  hyperbolic  trajectory  may  have  stable 
and  unstable  manifolds.  These  finite-length  manifolds  are  the 
material  lines  which  show  maximal  linear  stability /instability 
for  finite  times.  Hyperbolic  trajectories  and  manifolds  were 
shown  by  Duan  and  Wiggins  (1996)  and  Ide  et  al.  (2002)  to 
control  particle  transport  and  persist  in  time-dependent  flows. 
The  stable  and  unstable  manifolds  of  a  pair  of  hyperbolic 
fixed  points  were  found  by  Rogerson  (1999)  to  control  trans¬ 
port  in  and  out  of  a  cat’s  eye  region  for  a  barotropic  jet.  In 
genesis  cases,  the  time  of  hyperbolicity  is  short  enough  so 
that  multiple  intersections  of  manifolds  does  not  generally 
occur. 

Finite-time  Lagrangian  scalar-field  methods  were  devel¬ 
oped  (Haller  and  Poje,  1997;  Haller,  2001;  Haller  and  Yuan, 
2000;  Pierrehumbert,  1991)  as  an  efficient  way  of  directly 
locating  kinematic  boundaries  in  time-dependent  flows  with¬ 
out  computing  a  hyperbolic  trajectory.  These  studies  showed 


that  the  manifolds  were  the  material  lines  that  tended  to  have 
maximum  times  of  linear  instability.  The  manifolds  can  be 
located  as  maximal  ridges  in  Lagrangian  scalar  fields,  which 
measure  finite-time  separation  directly  from  particle  trajec¬ 
tories.  These  boundaries  will  hereafter  be  referred  to  as  La¬ 
grangian  coherent  structures  (LCSs).  LCSs  are  robust  under 
time- variation  of  fluid  velocities  (Shadden,  2006),  and  ve¬ 
locity  approximations  (Haller,  2002),  though  the  loss  of  hy¬ 
perbolicity  of  the  underlying  saddle  trajectory  may  alter  the 
invariance  of  the  LCS  (Branicki  and  Wiggins,  2010).  These 
methods  detect  LCSs,  which  form  finite-time  boundaries  and 
govern  the  mixing  and  transport  of  tracers  (Shadden  et  al., 
2005;  Haller  and  Yuan,  2000;  Haller,  2000)  and  have  proven 
useful  in  a  variety  of  geophysical  flows,  see  e.g.  Salman  et 
al.  (2008);  Tang  et  al.  (2009);  Koh  and  Legras  (2002);  Joseph 
and  Legras  (2002);  Cohen  and  Schultz  (2005). 

Many  powerful  diagnostic  tools  have  been  developed  by 
the  dynamical  systems  community  for  use  in  time-dependent 
flows.  However,  these  new  methods  are  not  well  known  to 
hurricane  forecasters  and  tropical  meteorologists. 

For  vortex  dominated  flows,  time- scale  finite-time  La¬ 
grangian  methods  were  introduced.  When  the  boundaries  are 
difficult  to  see  due  to  turbulence,  statistical  measures  of  tra¬ 
jectory  motion  give  insight  into  the  mixing  and  transport  pro¬ 
cesses. 

In  this  study,  Lagrangian  kinematic  aspects  of  the  pouch 
theory  are  explored  by  an  analysis  of  trajectory  motion  and 
scalar  fields  derived  from  trajectory  motion  using  methods 
from  dynamical  systems.  LCSs  are  located  as  the  maximal 
ridges  of  scalar  fields.  The  stirring,  which  is  modulated  by 
the  LCSs,  shows  the  pathways  for  the  entry  of  air  into  a  can¬ 
didate  disturbance  or  pre-storm. 

1.4  Summary 

The  outline  of  this  paper  is  as  follows.  In  Sect.  2,  we  provide 
an  introduction  to  Lagrangian  methods  for  time-dependent 
flows,  and  show  how  they  can  be  made  applicable  to  the  flows 
considered  during  the  PREDICT  experiment.  Section  3  pro¬ 
vides  numerical  details  of  the  Lagrangian  computations  as 
well  as  numerical  details  about  the  data  obtained  from  the 
PREDICT  experiment.  Sections  4  and  5  present  the  primary 
results  of  this  study.  In  Sect.  4,  we  consider  the  nature  of 
time-dependence  of  fluid  velocities  during  pre-genesis  condi¬ 
tions  and  contrast  Lagrangian  boundaries  to  Eulerian  stream¬ 
lines.  In  Sect.  5,  finite-time  Lagrangian  methods  and  La¬ 
grangian  statistical  methods  are  applied  to  a  developing  and 
a  non-developing  disturbance.  The  LCSs  are  shown  to  be  im¬ 
portant  for  controlling  the  import  of  dry  air  into  the  core  for  a 
non-developing  disturbance,  and  protecting  the  core  of  a  de¬ 
veloping  disturbance.  Conclusions  and  a  future  outlook  are 
provided  in  Sect.  6. 
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2  Overview  of  mathematical  methods 

The  methods  used  in  this  study  can  be  divided  into  measures 
of  recirculation  and  the  computation  of  kinematic  bound¬ 
aries.  For  both  types  of  methods,  we  consider  the  behavior  of 
particle  trajectories.  Let  xo  0^o(xo)  be  the  flow  map  asso¬ 
ciated  with  the  ordinary  differential  equation  for  fluid  veloc¬ 
ities 

x  =  m(x,  t)  (1) 

governing  the  location  of  particle  trajectories  x(xo,  to,t),  that 
is,  the  solution  of  Eq.  (1)  with  initial  condition  x(to)  =  xo- 
The  growth  of  perturbations  to  initial  conditions  can  be  ap¬ 
proximated  by  the  linearized  velocity  field, 

S  =  Vm§.  (2) 

2.1  OW 

Assuming  a  slowly  varying  flow,  the  growth  of  §  is  deter¬ 
mined  by  eigenvalues  of  Vm.  If  u  is  divergence  free,  the 
eigenvalues  of  Vm  are  opposite.  The  tendency  of  rotation 
of  an  air  parcel  is  determined  by  purely  imaginary  eigenval- 
ues  of  Vm,  k±  =  ±i  VOW,  when  OW=  \  (f: 2  -  -  Sj)  >  0. 
Sn  =  ux  —  vy  and  Ss  =  vx  +  uy  are  the  normal  and  parallel 
strain  rates,  f  =  vx  —  uy  is  vorticity,  and  i2  =  —  1.  In  this  sign 
convention,  positive  OW  indicates  that  vorticity  dominates 
strain.  Negative  OW  values  indicate  regions  of  high  strain, 
but  particles  in  regions  of  negative  OW  do  not  generally  ex¬ 
perience  persistent  stretching  even  in  a  steady  flow  (Lapeyre 
et  al.,  1999).  Hua  and  Klein  (2002)  provided  a  generalized 
OW  criterion  which  incorporates  time-dependent  terms.  An 
OW  criterion  incorporating  divergent  terms  is  studied  by 
Lukovich  and  Sheperd  (2005). 

2.1.1  Hyperbolic  fixed  point  and  manifolds 

Assume  that  a  velocity  field  «(x),xel2  is  steady.  A  fixed 
point  Xh  is  a  point  where  velocities  vanish,  u(xh)  =  0.  We 
are  concerned  only  with  fixed  points  having  hyperbolic  or 
saddle-type  stability.  Hyperbolicity  of  xu  may  be  defined  by 
Vm(x^)  having  eigenvalues  of  opposite  signs,  showing  a  di¬ 
rection  of  expansion  and  contraction.  The  stable  manifold 
is  the  solution  that  converges  to  the  hyperbolic  fixed  point 
along  the  eigendirection  of  the  positive  eigenvalue.  The  un¬ 
stable  manifold  diverges  from  the  hyperbolic  fixed  point  in 
the  eigendirection  of  the  negative  eigenvalue  or  converges 
backward  in  time.  The  stable  and  unstable  manifolds  are  im¬ 
portant  solutions  because  they  form  a  partition  of  the  flow, 
which  distinguish  them  from  generic  streamlines.  Particles 
near  the  stable  manifold  depart  and  travel  toward  the  unsta¬ 
ble  manifold  once  they  approach  the  fixed  point.  The  stable 
manifold  shows  maximal  repulsion  to  nearby  particles  while 
the  unstable  manifold  shows  maximal  attraction,  (Haller  and 
Yuan,  2000;  Ottino,  1990). 


2.2  Objectivity  of  reference  frames 

While  the  manifolds  of  a  hyperbolic  fixed  point  form  a  flow 
partition  that  is  conceptually  beneficial  for  describing  the 
pouch  boundary,  the  location  of  the  fixed  point  and  manifolds 
are  dependent  on  the  choice  of  translation  speed.  Moreover, 
the  rotating  flow  present  during  cyclogenesis  contributes  a 
non-Galilean  rotation  to  any  fixed  coordinate  system.  Thus 
flow  boundaries  in  the  co-moving  frame  may  be  significantly 
different  than  the  flow  boundaries  that  are  associated  with 
actual  particle  motion.  An  objective2  frame  of  reference  fol¬ 
lowing  particle  trajectories  avoids  many  of  the  approxima¬ 
tions  employed  previously  in  the  observational  and  model¬ 
ing  studies  testing  the  pouch  theory,  such  as  a  nearly  steady 
flow  in  a  co-moving  reference  frame.  In  addition,  objective 
measures  do  not  require  the  location  of  a  pouch  center  or  the 
computation  of  the  pouch  translation  speed. 

2.2.1  Time  variation,  hyperbolic  trajectories,  and  lobe 
dynamics 

Analogs  to  a  hyperbolic  fixed  point  and  its  invariant  mani¬ 
folds  may  still  exist  in  a  time-dependent  flow  and  have  the 
same  local  kinematic  properties  as  in  a  steady  flow.  A  distin¬ 
guished  hyperbolic  trajectory  (DHT)  is  a  particle  trajectory 
that  is  the  time-dependent  analog  of  a  hyperbolic  fixed  point 
and  maintains  linearized  hyperbolic  stability,  with  eigenval¬ 
ues  of  opposite  signs  at  each  time.  DHTs  are  often  marked 
by  persistent  hyperbolic  stagnation  points  at  each  fixed  time, 
but  the  persistence  of  a  hyperbolic  fixed  point  is  neither  a 
sufficient  or  necessary  condition  for  the  existence  of  a  hy¬ 
perbolic  trajectory.  (Branicki  and  Wiggins,  2010).  The  stable 
and  unstable  manifolds  of  the  DHT  are  the  particle  trajecto¬ 
ries  that  approach  the  DHT  forward/backward  in  time  (Ide  et 
al.,  2002).  The  manifolds  of  the  time-dependent  flow  are  not 
crossed  by  trajectories,  and  form  flow  boundaries  in  the  same 
manner  that  velocity  streamlines  do  for  steady  flows  (Malho- 
tra  and  Wiggins,  1999;  Coulliette  and  Wiggins,  2000).  A  key 
feature  of  manifolds  in  time-dependent  flows  is  that  the  en¬ 
tire  manifold  evolves  with  time.  Algorithms  by  Mancho  et 
al.  (2003)  and  Mancho  et  al.  (2004)  locate  manifolds,  but 
are  sensitive  to  the  persistence  of  stagnation  points  and  time- 
dependence  of  velocities.  Algorithms  by  Koh  and  Plumb 
(2000)  and  Koh  and  Legras  (2002)  are  not  sensitive  to  the 
persistence  of  stagnation  points.  Similar  to  the  steady  flow 
case,  the  stable  and  unstable  manifolds  intersect  at  the  hy¬ 
perbolic  trajectory. 


2 Objective  coordinate  transformations  are  of  the  form  x  i-> 
Q0)x  +  b(t ),  where  Q (t)  includes  rotations  of  the  coordinate  sys¬ 
tem.  These  transformations  are  more  general  than  Galilean  trans¬ 
formations  which  are  of  the  form  x  i->  Qx  +  ct  where  Q  is  a  proper 
orthogonal  matrix,  which  only  includes  translations.  For  a  thorough 
overview  of  the  definitions  and  choice  of  reference  frame  relating 
to  vortices,  see  Haller  (2005). 
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Stable  and  unstable  manifolds  are  allowed  to  intersect  at 
points  other  than  the  DHT.  If  additional  intersections  occur, 
contained  regions  formed  by  the  transverse  intersections  of 
the  manifolds  are  called  lobes.  Lobe  dynamics  describes  the 
transport  of  material  contained  within  the  lobes  (Malhotra 
and  Wiggins,  1999).  If  a  hyperbolic  trajectory  is  hyperbolic 
for  a  finite  time  interval,  the  manifolds  are  finite  length. 

2.3  Finite-time  Lagrangian  methods 

In  the  absence  of  persistent  stagnation  points  marking  a  DHT, 
finite-time  Lagrangian  methods  may  still  locate  manifolds. 
For  a  time-dependent  flow,  consider  two  nearby  trajectories 
straddling  a  manifold.  It  is  expected  that  these  trajectories 
experience  more  relative  separation  than  generic  pairs  of  tra¬ 
jectories.  Finite-time  Lagrangian  methods  measure  this  sep¬ 
aration  and  show  LCSs  which  correspond  to  manifolds  as  the 
regions  of  maximal  separation.  Though  the  DHT  is  not  com¬ 
puted  directly,  it  can  be  inferred  as  the  intersection  of  stable 
and  unstable  manifolds  located  near  a  persistent  hyperbolic 
fixed  point. 

2.3.1  Finite  time  Lyapunov  exponents 

Small  perturbations  in  the  initial  condition,  y0  =  xo  +  §0’ 
lead  to  a  perturbed  trajectory  y(t )  =  x(t )  +  $(0-  For  suffi¬ 
ciently  small  |§ol»  the  perturbation  §  (t)  can  be  approximated 
through  the  Jacobian  of  the  flow  map  as 

SO)  =dxo<pl,0  (*o)£o>  (3) 

which  satisfies  the  variational  equation 

$  =  Vii(jr(f),f)$.  (4) 

Consider  a  time  range  [to,  to  +  8 1]  with  fixed  integration  time 
8t  which  is  finite  but  not  necessarily  small.  The  growth  of 
£(7)  during  at  time  t  +8t  \ s  approximated  by 

IISO)ll2  =  (So,  {dXo^0+s\x0))*(dxJrt°0+Sr(x0))^).  (5) 

where  *  denotes  a  matrix  transpose.  The  Cauchy-Green  de¬ 

formation  tensor  defined  as 

A(x0,  t0,  t0 +  8t)  =  (dxo(/)l00+s\xo))* (dX0<p‘t°+St (x0)). 

Maximal  growth  of  the  perturbation  occurs  when  §(0)  is 
aligned  with  the  eigenvector  corresponding  to  the  larger 
eigenvalue,  A.max(A),  of  A, 

max  US  (Oil  =  \/<So.  Vax(  A(x0,  f<).  *o  +  <$0)S  (0))  (6) 

So 

=  y/ ^max(A(Xo,  ^0^0  +  ^0)  II §0  II  •  (7) 

The  objective  quantity 

1 

a(xo,to,to  +  &t)  =  — —  lnkmax(A(A:o^o^o  +  50) 

2\bt\ 

is  the  FTLE  (considered  as  a  function  of  xo)  at  the  initial 
time  to. 


2.3.2  Lagrangian  scalar  fields 

Finite-time  Lagrangian  scalar-field  methods  compute  a  field 
of  Lagrangian  values  on  a  grid  of  initial  conditions.  La¬ 
grangian  values  are  assigned  the  initial  position  xo-  Note  that 
for  finite-times,  the  parameter  8t  =  t  —  to  is  the  chosen  inte¬ 
gration  time.  The  scalar  field  of  FTLEs  is  computed  by  eval¬ 
uating  the  FTLE  on  a  grid  of  initial  conditions.  Repelling  and 
attracting  LCSs  are  maximal  ridges  of  the  forward  (8t  >  0) 
and  backward  ( 8t  <  0)  FTLE-fields,  which  play  the  role  of 
stable  and  unstable  manifolds  over  the  finite  integration  time 
(Haller,  2000;  Shadden  et  al.,  2005;  Shadden,  2006;  Haller, 
2002). 

3  Data  from  the  PREDICT  Experiment 
3.1  Analysis  data 

The  data  used  for  this  study  from  the  PREDICT  experiment 
is  taken  from  the  ECMWF  model,  which  outputs  velocity 
and  moisture  data  on  an  evenly  spaced  planar  grid  with  a 
spatial  resolution  of  0.25  degrees  and  temporal  output  every 
12  h.  Vertical  levels  for  velocity  data  output  are  on  constant 
pressure  levels.  In  this  study,  the  ECMWF  model  output  and 
forecast  data  for  the  case  of  two  tropical  disturbances  from 
28  August  to  17  September  is  used,  during  the  times  that 
Gaston  and  Karl  were  tracked.  The  disturbance  that  became 
Gaston  had  a  visible  pouch  on  28  August,  and  was  designated 
a  tropical  storm  by  National  Hurricane  Center  forecasters  on 
30  August  and  was  later  downgraded  on  2  September  to  an 
area  of  disturbed  weather.  The  Ex-Gaston  was  trackable  until 
1 1  September,  but  never  redeveloped.  The  disturbance  that 
eventually  became  Karl  had  a  well-defined  surface  low  along 
the  wave  trough  and  associated  pouch  on  9  September,  and 
lasted  until  17  September.  Further  details  on  the  evolution  of 
these  storms  is  given  in  Montgomery  et  al.  (2012),  while  the 
Lagrangian  behaviour  of  these  storms  is  described  in  Sect.  4. 

The  temporal  output  of  forecast  data  is  every  12  h,  with  a 
120h  forecast  run  each  day.  The  use  of  Lagrangian  methods 
to  locate  structures  in  the  forecast  data  would  reveal  LCSs 
that  would  exist  under  the  given  forecast.  The  forecast  data 
provides  a  way  to  compute  repelling  LCSs  in  real  time,  but 
will  be  addressed  in  a  future  study.  For  this  study,  we  use 
only  the  available  analysis  output  data  at  Oh  and  12 h  each 
day. 

This  study  focuses  on  the  kinematic  aspects  of  genesis, 
which  when  combined  with  the  dynamic  and  thermodynamic 
aspects,  see  Smith  and  Montgomery  (2012),  will  help  pro¬ 
vide  a  more  complete  picture  of  genesis  in  the  real  atmo¬ 
sphere.  For  the  kinematic  boundaries,  we  focus  on  the  model 
analysis  data  which  incorporates  dropsonde  data,  but  we  do 
not  use  dropsonde  data  directly  since  the  Lagrangian  compu¬ 
tations  require  trajectory  integrations  which  are  far  simpler 
using  the  model  output  grids. 
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3.2  Trajectory  computations 

All  of  the  Lagrangian  fields  are  computed  daily  in  a  24  de¬ 
gree  longitude  by  1 8  degree  latitude  domain  centered  at  the 
pouch  center  by  seeding  particles  at  0.125  degree  intervals 
and  computing  particle  trajectories  both  forward  and  back¬ 
ward  in  time.  The  time  interval  for  model  data  is  12  h,  but 
we  use  a  15  minute  intermediate  time  step  on  a  fourth  or¬ 
der  Runge-Kutta  solver  to  better  resolve  trajectory  curvature. 
Decreasing  the  intermediate  time-step  further  showed  no  vis¬ 
ible  effects  on  Lagrangian  fields.  Due  to  the  coarse  temporal 
output  and  intrinsic  curvature  of  the  flow,  we  adopt  a  simple 
feature  tracking  scheme  for  the  temporal  interpolation  of  ve¬ 
locities  at  intermediate  times.  The  scheme  aligns  particles  at 
intermediate  times  to  their  feature-relative  position  at  known 
model  output  times  by  advancing  the  entire  feature  at  the 
mean-flow  speed.  For  a  particle  x(t)  and  t  an  intermediate 
time,  the  velocity  of  the  particle  at  the  known  output  time  4 
is  u(x(tk))  =  u(x(t)  +  cp(t  —  tk))  where  cp  is  the  computed 
translation  speed.  Linear  interpolation  in  time  is  then  used  to 
obtain  fluid  velocities  at  the  intermediate  output  times.  Lin¬ 
ear  spatial  interpolation  is  used  to  compute  velocities  after 
the  feature  alignment  is  done.  This  process  improves  on  sim¬ 
ple  linear  or  bi-cubic  interpolation  since  curved  trajectories 
can  be  tracked  through  a  translating  vortex. 

In  order  to  balance  the  constraints  of  resolution,  availabil¬ 
ity  of  information,  and  computational  expense,  we  have  cho¬ 
sen  48  h  as  the  integration  time.  This  time  is  long  enough 
to  resolve  individual  LCSs  while  maintaining  coherence  of 
LCS’s  across  varying  initial  times.  Shorter  integration  times 
did  not  resolve  LCS’s  as  well,  while  much  longer  integra¬ 
tion  times  resolved  LCSs  but  were  computationally  expen¬ 
sive  and  sometimes  yielded  spurious  structures  associated 
with  the  domain  boundary.  The  occurrence  of  these  spuri¬ 
ous  structures  was  reduced  by  allowing  particles  that  con¬ 
tacted  the  southern  domain  boundary  to  travel  zonally  at  the 
same  speed  as  particles  at  the  equator  while  further  pole- 
ward  motion  is  suppressed.  Since  they  are  based  on  parti¬ 
cle  trajectories,  Lagrangian  fields  require  48  h  of  velocities, 
but  only  require  the  computed  pouch  location  as  a  point  of 
reference  for  visualization  at  the  initial  time  of  integration. 
As  explained  in  Sect.  3.3,  we  assume  that  the  flow  is  layer- 
wise  two-dimensional,  and  have  chosen  500  hPa,  600  hPa, 
700  hPa,  850  hPa,  and  925  hPa  as  vertical  levels  for  all  com¬ 
putations. 

Vorticity  and  OW  values  require  the  computation  of  ve¬ 
locity  gradients,  which  are  computed  through  second-order 
finite  differences. 

3.3  Simplifications  used  for  this  study 

During  the  PREDICT  experiment,  for  each  model  output 
time,  a  pouch  translation  speed  was  computed  to  view 
the  co-moving  reference  frame.  Particle  trajectories  fol¬ 
low  co-moving  streamlines  if  the  flow  is  steady  in  the  co¬ 


moving  frame  and  layer-wise  two-dimensional.  Finite-time 
Lagrangian  methods  require  only  the  computation  of  particle 
trajectories,  so  the  assumptions  of  a  steady  flow  is  unnec¬ 
essary.  Lagrangian  structures  are  found  independent  of  the 
co-moving  frame  of  reference  associated  with  a  computed 
translation  speed.  However,  the  translation  speed  is  still  used 
to  compute  velocities  at  times  between  model  output  times 
as  described  above. 

Lagrangian  methods  have  been  developed  that  are  suitable 
for  three-dimensional  flows,  (Haller,  2005,  2011;  Ruther¬ 
ford  and  Dangelmayr,  2010).  However,  the  interpretation 
of  boundaries  is  easier  also  in  two  dimensions  since  the 
manifolds  are  one-dimensional.  Justification  for  computing 
two-dimensional  structures  on  horizontal  layers  in  a  three- 
dimensional  flow  is  provided  by  Rutherford  et  al.  (2012), 
which  showed  that  the  two-dimensional  projection  of  LCSs 
computed  from  three-dimensional  trajectories  were  aligned 
across  vertical  levels  up  to  7  km  in  the  vicinity  of  rotating 
deep  convective  structures;  the  flow  geometry  was  nearly 
two-dimensional  even  though  trajectory  paths  were  not.  A 
similar  alignment  of  LCSs  across  vertical  levels  is  shown  in 
this  study. 

Despite  the  availability  of  3-D  Lagrangian  methods,  these 
methods  can  only  be  employed  if  3-D  trajectory  locations 
are  accurate.  We  use  the  simplification  of  a  2-D  velocity  field 
due  to  the  higher  reliability  of  2-D  motion  from  the  avail¬ 
able  global  model  data.  Trajectories  from  2-D  velocities  are 
suitable  for  analysing  large-scale  motions  associated  with  the 
primary  pouch  boundary  since  the  tropical  atmosphere  is  typ¬ 
ically  stratified  except  in  convective  regions  (Holton,  2004), 
which  typically  occur  inside  the  pouch  boundary.  A  layer- 
wise  two-dimensional  flow  is  further  justified  since  convec¬ 
tive  features  are  not  well  resolved  by  global  models,  thus 
horizontal  convergence/divergence  is  small.  In  convective  re¬ 
gions,  ascent  times  are  on  the  order  of  1  h,  which  is  much 
smaller  than  the  gap  between  model  output  times.  The  hor¬ 
izontal  width  scale  of  strong  convective  features  is  approxi¬ 
mately  10  km  which  is  smaller  than  the  highest  global  model 
output  resolution  of  25  km.  (Houze  Jr.  et  al.,  2009;  Bell  and 
Montgomery,  2010;  Reasor,  2005) 

Since  we  use  a  trajectory-based  approach,  the  “vertical 
escape”  of  trajectories  off  of  the  horizontal  slices  is  a  con¬ 
cern  of  the  simplification.  While  trajectories  may  escape,  the 
persistence  of  coherent  structures  across  vertical  levels  sug¬ 
gests  that  the  flow  geometry  remains  approximately  planar 
even  though  trajectory  motions  may  not.  Approximate  ver¬ 
tical  alignment  of  LCSs  is  required  for  the  LCSs  to  form  a 
two-dimensional  manifold  in  the  three-dimensional  flow. 


4  An  analysis  of  time-dependence  in  tropical 
disturbances 

In  this  section,  we  examine  the  kinematic  structures  which 
form  finite-time  boundaries  and  control  particle  transport 
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between  the  environment  and  the  pouch.  Examples  of  these 
boundaries  will  be  shown  for  ex-Gaston  while  it  was  clas¬ 
sified  as  a  tropical  disturbance.  We  find  that  while  the  La¬ 
grangian  boundaries  are  related  to  the  Eulerian  manifolds, 
time-dependence  and  finite-time  hyperbolicity  cause  impor¬ 
tant  differences. 

4.1  Reference  frames 

In  the  co-moving  frame,  a  clearer  view  of  the  geometry  asso¬ 
ciated  with  recirculation  can  be  seen  by  a  saddle  point  often 
appearing  outside  of  the  pouch,  and  stable  and  unstable  man¬ 
ifolds  forming  the  pouch  boundary,  Fig.  1.  The  unstable  man¬ 
ifold  encloses  the  pouch,  while  the  region  between  the  stable 
and  unstable  manifolds  is  a  pathway  for  environmental  air  to 
reach  the  circulation  center.  In  the  case  that  the  stable  and 
unstable  manifolds  of  a  single  fixed  point  coincide,  a  tear¬ 
drop  shaped  boundary  is  formed  and  blocks  any  transport 
into  or  out  of  the  pouch.  If  the  manifolds  of  a  pair  of  fixed 
points  coincide,  a  cat’s  eye  boundary  is  formed.  Figure  2b 
shows  the  co-moving  streamlines  and  a  cat’s  eye  boundary, 
along  with  trajectory  locations  which  enter  a  circle  of  3  de¬ 
gree  radius  around  the  pouch  center  (yellow)3  within  48  h  at 
the  same  time  that  the  streamlines  are  shown.  Some  parti¬ 
cles  that  enter  the  pouch  from  the  northeast  not  only  cross 
Eulerian  streamlines,  Fig.  2a,  but  also  cross  the  co-moving 
streamlines,  Fig.  2b.  Since  the  3  degree  circle  is  completely 
contained  within  the  co-moving  cat’s  eye,  the  trajectories 
shown  do  not  represent  all  of  the  air  that  enters  the  pouch 
from  outside  the  cat’s  eye.  The  co-moving  frame  does  not 
account  for  all  2-D  particle  motion  for  two  reasons.  The  first 
reason  is  that  the  co-moving  frame  is  based  on  a  translation 
speed  that  is  assumed  constant  in  both  space  and  time.  The 
second  reason  is  that  velocities  in  the  co-moving  frame  are 
still  time-dependent. 

4.2  Sensitivity  to  changes  in  translation  speed 

Recirculation  is  the  primary  characteristic  of  a  vortex  which 
enables  repeated  moistening  and  vorticity  amplification,  and 
is  generally  not  sensitive  to  small  variations  in  the  computed 
pouch  translation  speed.  However,  the  flow  boundaries  out¬ 
side  the  recirculation  region  directly  impact  which  environ¬ 
mental  air  masses  interact  with  the  vortex  and  may  be  sen¬ 
sitive  to  variations  in  the  translation  speed.  A  slight  change 
in  translation  speed  may  change  the  shape  and  location  of 
manifolds.  Time- variations  and  spatial  variations  between 
the  inner-core  circulation  and  the  hyperbolic  fixed  point  lo¬ 
cation  of  the  pouch  translation  speed  may  alter  the  streamline 

3  A  yellow  circle  in  a  figure  represents  a  3  degree  radius  circle 
about  the  pouch  center,  which  was  computed  as  the  location  of  the 
intersection  of  the  wave  critical  layer  and  wave  trough  axis.  PGI38L 
and  PGI44L  refer  to  Gaston  and  Karl  respectively.  Longitude  is 
given  on  a  0-360°  convention.  In  this  convention,  the  longitude  of 
St.  Croix  is  295.25°. 


Fig.  1.  Schematic  of  the  hyperbolic  fixed  point  and  manifolds  near 
a  convergent  vortex  in  a  layer- wise  two-dimensional  flow. 


structure  in  the  co-moving  frame.  The  time  variation  of  the 
computed  pouch  translation  speed  over  a  12  hour  time  inter¬ 
val  may  exceed  2  ms-1.  Spatial  variations  in  velocities  be¬ 
tween  a  pouch  center  and  pouch  boundary  also  often  exceed 
2ms-1. 

The  effect  of  a  spatially  uniform  2ms-1  change  in  transla¬ 
tion  speed  to  streamline  locations  is  shown  in  Fig.  2c.  There 
are  saddle  points  to  the  northeast  and  west  of  the  distur¬ 
bance,  and  the  northeast  saddle  point  and  associated  stream¬ 
lines  constrains  the  movement  of  the  particles  to  the  north. 
The  co-moving  streamlines  indicate  that  the  particles  to  the 
north  will  move  south  toward  the  storm,  but  will  be  taken 
to  the  east  as  they  move  southward,  and  will  not  enter  the 
storm.  Streamlines  computed  with  the  different  translation 
speed  show  a  much  larger  pouch  boundary  and  the  manifolds 
of  the  western  stagnation  point  do  not  form  a  complete  en¬ 
closure  around  the  pouch.  The  streamlines  indicate  that  the 
particles  from  the  north  will  now  be  recirculated  within  the 
pouch.  In  this  example,  the  particle  motion  is  more  accu¬ 
rately  described  by  the  streamlines  computed  with  the  ad¬ 
justed  translation  speed.  In  other  cases,  the  particles’  final 
locations  are  mistaken  by  the  streamlines  computed  from  ei¬ 
ther  translation  speed,  where  in  fact,  many  particles  outside 
the  boundary  enter  the  pouch.  Since  the  particles  originate 
from  a  region  with  low  relative  humidity,  the  streamlines  un¬ 
derestimate  the  amount  of  dry  air  entering  the  pouch.  The 
differences  in  translation  speed  tend  to  have  less  effect  on 
the  geometry  of  strong  vortices  (Riemer  and  Montgomery, 
2011).  Since  the  Lagrangian  frame  does  not  require  a  transla¬ 
tion  speed,  it  is  not  sensitive  to  variations  in  translation  speed. 

4.3  LCSs 

LCSs  are  the  finite-time  approximations  to  stable  and  un¬ 
stable  manifolds  of  a  hyperbolic  trajectory.  In  Fig.  2d,  the 
forward  (red)  and  backward  time  (blue)  FTLE  ridges  for  ex- 
Gaston  are  shown  at  700  hPa.  There  is  a  prominent  attract¬ 
ing  LCS  (blue)  and  repelling  LCS  (red)  which  are  in  nearly 
the  same  location  as  the  unstable  and  stable  manifolds  of  the 
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Fig.  2.  Streamlines  of  Earth-relative  velocities  (a),  streamlines  of  co-moving  velocities  (b),  and  streamlines  of  comoving  velocities  (c)  with  a 
+2  m  s_1  variation  in  translation  speed  shown  on  2  September  for  ex-Gaston  and  overlaid  on  the  relative  humidity  field  (in  %)  (a-c)  Particle 
locations  beginning  in  the  pouch  on  4  September  are  shown  as  green  dots  in  (b).  Co-moving  stable  and  unstable  manifolds  are  shown  as  red 
and  blue  lines,  respectively,  in  (b)  and  (c).  Attracting  (blue)  and  repelling  (red)  LCSs  are  represented  as  level  contours  of  the  FTLE  field  at 
80  %  of  the  maximal  FTLE  values,  and  are  overlaid  on  the  relative  humidity  field  (d).  The  attracting  LCSs  are  overlaid  on  the  duration  time 
of  hyperbolic  OW  (hours)  and  duration  time  of  elliptic  OW  (hours)  in  (e)  and  (f)  respectively. 
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eastern  stagnation  point  of  the  co-moving  frame.  Unlike  the 
streamlines,  the  LCSs  have  finite  length,  and  do  not  connect 
with  the  western  stagnation  point.  The  LCSs  are  not  crossed 
by  particle  trajectories.  To  reach  the  pouch  center,  trajecto¬ 
ries  travel  around  the  LCS,  which  allows  the  set  of  trajec¬ 
tories  from  the  northeast  to  enter  the  core.  The  LCS  wraps 
around  the  core,  allowing  a  recirculation  region  inside  the 
LCS.  In  contrast  to  the  eastern  stagnation  point,  the  mani¬ 
folds  of  the  western  stagnation  point  are  not  represented  by 
persistent  LCSs.  There  are  some  very  short  LCS  segments 
where  the  manifolds  are  located,  but  they  do  not  connect  with 
the  manifolds  of  the  eastern  stagnation  point  to  form  an  im¬ 
penetrable  boundary. 

Differences  in  the  LCSs  and  streamlines  away  from  the 
hyperbolic  fixed  point  are  due  to  the  time-variation  in  ve¬ 
locities.  However,  if  the  velocity  near  a  stagnation  point  is 
slow,  the  co-moving  manifolds  may  be  computed  while  the 
LCS  associated  with  the  manifold  does  not  exist,  since  there 
is  not  sufficient  particle  separation.  The  nonexistence  of  an 
LCS  where  a  stable  manifold  resides  provides  one  explana¬ 
tion  to  how  particles  appear  to  pass  through  the  co-moving 
manifold. 

FTLE  values  along  repelling  LCSs  during  periods  of  con¬ 
vergence  are  often  smaller  than  the  FTLE  values  along  the 
attracting  LCS.  While  LCSs  with  high  FTLE  values  tend  to 
be  well  resolved  and  persistent  across  varying  initial  time,  re¬ 
pelling  LCSs  with  smaller  FTLE  values  do  not  persist  across 
varying  initial  times.  Small,  isolated  ridges  with  low  FTLE 
values  that  do  not  persist  across  varying  initial  times  are  arti¬ 
facts  of  the  Eulerian  flow  that  would  disappear  under  a  longer 
integration  time,  but  are  still  visible  due  to  finite  integration 
times  where  the  velocity  at  a  particular  output  time  still  has 
some  influence  over  trajectory  locations. 

The  relation  between  the  FTLE  values  and  the  time  that 
particles  are  hyperbolic  can  be  seen  in  Fig.  2e,  where  the 
attracting  LCS  is  overlaid  on  the  field  showing  the  number 
of  hours  in  a  72  h  backward  integration  that  particles  satisfy 
the  hyperbolic  condition  of  the  OW  criterion.  The  strongest 
ridges  are  those  where  particles  remain  hyperbolic  for  the 
longest  times,  (Haller  and  Yuan,  2000).  For  this  reason,  we 
discard  LCSs  with  low  FTLE  values  since  they  tend  to  not 
be  attracting  or  repelling  to  nearby  particles  for  long  periods 
of  time,  and  would  not  be  expected  to  have  any  meaningful 
influence  as  flow  boundaries. 

If  attracting  and  repelling  LCSs  of  the  eastern  and  west¬ 
ern  hyperbolic  trajectories,  respectively,  are  resolved,  and  the 
LCSs  have  multiple  intersections  to  the  north  of  the  pouch 
center  near  the  moist  and  dry  air  boundary,  transport  across 
the  Eulerian  manifold  is  still  possible  through  lobe  dynamics. 
Enclosed  lobes  formed  from  intersecting  segments  of  LCSs 
allow  material  to  enter  the  circulation  center  through  the  de¬ 
formation  of  the  LCSs  even  when  the  Eulerian  manifolds  are 
persistent,  though  the  air  does  not  cross  the  LCSs.  For  the 
case  of  Gaston,  lobe  dynamics  did  not  occur  due  to  the  fi¬ 
nite  hyperbolicity  time  and  manifolds  with  short  lengths  that 


did  not  have  multiple  intersections.  For  Karl,  lobe  dynamics 
occurs  well  after  genesis. 

The  foregoing  considerations  offer  two  explanations  as  to 
how  the  particles  pass  through  the  co-moving  manifolds  to 
bring  dry  air  to  the  circulation  center. 

The  relation  of  LCSs  to  physical  quantities  is  shown  in 
Fig.  3.  Attracting  LCSs  define  a  sharp  contrast  between  sets 
of  particles  based  on  the  lengths  of  trajectory  paths,  Fig.  3a. 
In  this  case,  the  particles  north  of  the  LCS  take  much  longer 
paths  than  those  to  the  south  of  the  LCS.  Attracting  LCSs 
are  located  along  the  gradients  of  physical  quantities  that  are 
close  to  being  materially  conserved,  such  as  ozone,  Fig.  3b, 
and  0e,  Fig.  3c.  Particles  with  the  highest  potential  vorticity 
and  those  satisfying  the  elliptic  OW  criterion  tend  to  aggre¬ 
gate  at  the  interior  of  the  attracting  LCS,  Fig.  3d,  e.  In  con¬ 
trast  to  the  Eulerian  unstable  manifold,  the  attracting  LCS 
forms  both  the  boundary  that  protects  the  vortex  from  lateral 
intrusion,  and  the  interior  of  its  curve  marks  the  focal  point 
of  vorticity  aggregation. 

5  Analysis  of  a  non-developing  and  developing 
disturbance 

We  now  apply  the  Lagrangian  methods  to  the  cases  of  Gaston 
and  Karl  to  show  what  additional  insight  these  methods  can 
provide  in  describing  the  development  or  non-development 
of  a  tropical  depression. 

5.1  Evolution  of  ex-Gaston 

5.1.1  Summary 

To  provide  insight  as  to  why  ex-Gaston  did  not  re-develop, 
we  look  at  the  interaction  of  the  storm  with  its  environment  at 
700  hPa  and  above  which  caused  the  pouch  to  become  shal¬ 
lower  despite  a  favourable  environment  at  lower  levels.  A 
time  series  of  pouch  averaged  metrics  showing  the  evolution 
of  Gaston  is  shown  in  Fig.  4.  Gaston  was  declared  a  tropi¬ 
cal  storm  on  1  September  but  was  downgraded  to  a  tropical 
depression  on  2  September.  On  1  September,  Gaston  resided 
in  a  relatively  moist  environment  with  high  moisture  at  lower 
levels  but  surrounded  by  dry  air  at  upper  levels.  The  influence 
of  the  dry  air  can  be  seen  by  deep  shear  and  pouch  shear  that 
is  above  12 ms-1  prior  to  1  September,  and  remains  above 
8 ms-1  on  and  after  1  September.  Moisture  values  peaked 
on  1  September,  and  began  to  decline  at  700  hPa  and  above, 
as  substantial  dry  air  had  penetrated  the  pouch  at  mid  and 
upper  levels,  lowering  the  pouch-averaged  relative  humidity 
by  approximately  20%  by  5  September.  Before  5  Septem¬ 
ber,  this  intrusion  was  not  noticeable  at  700  hPa  in  the  rela¬ 
tive  humidity  analysis  data,  and  lower  levels  actually  showed 
moistening  (Montgomery  et  al.,  2012).  The  interaction  of  dry 
air  with  the  pouch  suggests  that  the  pouch  was  subjected  to 
vertical  wind  shear  (Davis  and  Ahijevych,  2011).  Smith  and 
Montgomery  (2012)  show  that  the  dry  air  entering  the  pouch 
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Fig.  3.  The  length  of  backward  trajectory  paths  in  degrees  (a),  O3  mass  fraction  (kg  kg-1)  (b),  6e  (K)  (c),  potential  vorticity  (s~l)  (d),  and 
OW  (s-2)  (e)  are  shown  for  Gaston  at  00:00  UTC  on  2  September  at  700  hPa  with  attracting  (blue)  and  repelling  (red)  LCSs  are  represented 
as  level  contours  of  the  FTLE  field  at  80  %  of  the  maximal  FTLE  values. 
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decreased  0e  in  the  mid  to  lower  troposphere  but  did  not 
produce  widespread  downdrafts  strong  enough  to  decrease 
the  boundary  layer  0e  on  the  pouch  scale.  On  1  September, 
the  pouch  has  high  relative  humidity  at  all  levels,  but  is  sur¬ 
rounded  by  dry  air.  As  the  humidity  decreased,  the  pouch- 
averaged  OW  and  vorticity,  Fig.  4c,  e,  also  began  to  decline. 

A  time  series  of  pouch- averaged  relative  humidity,  vor¬ 
ticity  and  OW  are  shown  in  Fig.  4.  Relative  humidity  val¬ 
ues  above  80  %  until  4  September,  Fig.  4a,  implies  that  dry 
air  entering  the  pouch  occurs  at  700  hPa  and  above  and  is 
compensated  by  convection  and  the  associated  vertical  trans¬ 
port  of  moisture  from  below.  This  view  is  supported  also  by 
the  study  of  the  thermodynamics  of  ex-Gaston  by  Smith  and 
Montgomery  (2012). 

The  Lagrangian  coherent  structures  and  trajectory  analysis 
will  demonstrate  that  the  intrusion  of  dry  air  was  occurring 
laterally  in  the  period  from  1  to  5  September,  but  the  air  was 
presumably  moistened  from  convection.  However,  the  pouch 
structure  was  significantly  degraded.  5  September  marked 
the  point  where  the  intrusion  of  dry  air  and  escape  of  mois¬ 
ture  through  ventilation  could  not  be  overcome  by  convective 
moistening  from  below.  By  this  point,  the  pouch  had  already 
lost  most  of  its  vertical  structure,  as  there  was  no  upper  level 
circulation  remaining  (Evans  et  al.,  201 1).  Time  series  show¬ 
ing  the  evolution  of  forward  and  backward  time  Lagrangian 
quantities  are  shown  in  Fig.  5. 

5.1.2  Analysis  of  the  flow  field 

The  co-moving  manifolds  are  shown  for  Gaston  from  1  to  6 
September  at  700  hPa  in  Fig.  6.  A  persistent  hyperbolic  stag¬ 
nation  point  can  be  seen  to  the  east  of  the  pouch  at  700  hPa, 
though  the  manifolds  do  not  completely  protect  the  pouch 
center.  On  2  September,  an  additional  stagnation  point  ap¬ 
pears  to  the  west  of  the  pouch,  and  a  cat’s  eye  configuration 
remains  in  place  through  6  September,  though  it  opens  tem¬ 
porarily  on  3  September. 

The  evolution  of  the  Lagrangian  fields  at  700  hPa  is  shown 
for  1  September  through  6  September  in  Fig.  7,  where  attract¬ 
ing  LCSs  (blue)  and  repelling  LCSs  (red)  are  overlaid  on  the 
relative  humidity  field.  There  are  LCSs  present  that  are  re¬ 
lated  to  the  co-moving  manifolds  of  the  eastern  hyperbolic 
stagnation  points.  These  structures  are  visible  on  1  Septem¬ 
ber,  Fig.  7a,  where  an  attracting  LCS,  shown  in  blue,  is  lo¬ 
cated  in  a  similar  location  as  the  unstable  manifold  (red)  of 
the  eastern  stagnation  point  of  the  cat’s  eye.  This  LCS  has 
a  tail  southeast  of  the  pouch  and  wraps  cyclonically  around 
the  north  side  of  the  pouch,  and  forms  the  northern  boundary. 
Near  the  location  of  the  eastern  stagnation  point,  a  repelling 
LCS  intersects  the  attracting  LCS,  and  the  intersection  marks 
a  hyperbolic  trajectory  that  is  persistent  until  5  September 
when  the  manifolds  no  longer  intersect.  A  repelling  LCS  is 
persistent  and  splits  the  large  region  of  high  moisture  on  the 
south  side  of  the  pouch,  eliminating  a  large  portion  of  that 


moisture  from  contributing  to  development  through  import¬ 
ing  moist  air. 

A  notable  difference  between  the  LCSs  and  the  Eulerian 
manifolds  is  that  the  LCSs  are  of  finite  length  because  of  the 
finite-time  hyperbolicity  of  the  hyperbolic  trajectory.  Since 
they  are  of  finite  length,  they  do  not  form  a  complete  enclo¬ 
sure  around  the  pouch  of  ex-Gaston. 

A  repelling  LCS  related  to  the  western  stagnation  point 
does  not  intersect  the  attracting  LCS  at  the  northern  edge  of 
the  pouch  to  form  a  cat’s  eye  like  boundary  from  air  masses 
to  the  north,  Fig.  7a-f.  Since  this  boundary  is  not  closed  and 
the  repelling  LCS  is  not  persistent,  the  attracting  LCS  con¬ 
trols  the  interaction  of  the  pouch  with  dry  air  to  the  north. 
The  attracting  LCS  forms  most  of  the  boundary  to  the  north, 
and  resides  approximately  3  degrees  from  the  pouch  center. 
The  attracting  LCS  forces  air  that  is  east  of  the  structure  to 
wrap  around  the  pouch,  and  the  air  does  not  enter  for  at  least 
two  days. 

Vertical  wind  shear  pushes  particles  from  a  broad  region  of 
low  relative  humidity  from  the  north  toward  the  pouch  (Davis 
and  Ahijevych,  2011),  and  forces  the  LCS  southward.  This 
dry  air  does  not  directly  enter  the  pouch  as  it  is  blocked  by 
the  attracting  LCS,  which  forces  the  dry  air  mass  to  wrap 
around  the  west  side  of  the  pouch  before  it  is  entrained.  The 
orientation  of  the  repelling  LCS  indicates  that  the  moist  air 
to  the  south  of  the  pouch  will  not  be  drawn  inward. 

By  2  September,  the  attracting  LCS  splits  east  of  the  pouch 
center,  begins  to  effectively  “unzip”,  Fig.  7d,  and  part  of  the 
LCS  begins  to  travel  cyclonically  around  the  center.  This 
branch  rejoins  the  stationary  branch  on  4  September  after 
making  one  revolution  around  the  pouch,  Fig.  7d,  further 
constricting  the  moisture  source  from  the  ITCZ,  Fig.  7d.  The 
western  repelling  LCS  rotates  with  the  cyclonic  flow  toward 
the  pouch  center,  and  does  not  prohibit  dry  air  from  enter¬ 
ing  since  it  is  not  constrained  by  multiple  intersection  points. 
As  it  enters  the  pouch  center,  it  disappears,  and  is  replaced 
by  a  new  repelling  LCS  to  the  north  of  the  pouch  center 
on  2  September,  Fig.  7c.  The  replacement  of  the  repelling 
LCS  allows  a  region  of  dry  air  between  the  repelling  LCS 
branches  to  cross  where  the  co-moving  streamline  resides 
from  2  September  to  4  September. 

A  pathway  for  material  exiting  the  pouch  also  emerges  due 
to  a  change  in  the  flow  field  and  is  marked  by  the  presence 
of  a  repelling  LCS,  which  is  entrained  into  the  pouch  cen¬ 
ter  at  700  hPa  and  above  on  4  September.  At  the  same  time, 
the  attracting  LCS  no  longer  encircles  the  pouch  center.  The 
aggregate  contribution  of  the  repelling  LCS  can  be  seen  by 
the  time  series  of  forward  time  FTLE  values  in  the  Figure  5, 
which  shows  higher  FTLE  values  at  500  and  700  hPa  begin¬ 
ning  on  2  September  at  the  time  where  significant  mixing 
was  occurring  in  the  pouch.  Repelling  LCSs  at  the  pouch 
center  normally  have  smaller  FTLE  values  but  this  LCS  has 
FTLE  values  similar  to  those  of  attracting  LCSs,  indicating  a 
shift  from  convergence  to  divergence  within  the  pouch.  After 
the  repelling  LCS  appears  and  the  ITCZ  moisture  source  is 
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Fig.  4.  Pouch  averages  of  Eulerian  quantities  for  the  2  degree  circle  surrounding  the  pouch  are  shown  for  ex-Gaston  (left)  and  pre-Karl  (right) 
for  the  relative  humidity  (%)  (a,  b),  £  (s_1)  (c,  d),  and  Okubo-Weiss  parameter  (s-2)  (e,  f).  The  different  curves  represent  vertical  levels  of 
500  hPa,  600  hPa,  700  hPa,  850  hPa,  and  925  hPa. 
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forward  time  FTLEs  forward  time  FTLEs 


Fig.  5.  Pouch  averages  of  backward-time  (a,  b)  and  forward  time  (c,  d)  FTLE  fields  for  a  2  degree  circle  surrounding  the  pouch  are  shown 
for  Gaston  (a,  c)  and  Karl  (b,  d).  The  different  curves  represent  vertical  levels  of  500  hPa,  600  hPa,  700  hPa,  850  hPa,  and  925  hPa. 


removed,  air  from  the  pouch  is  allowed  to  leave  the  pouch  to 
the  west. 

Trajectory  analysis 

The  relation  between  LCSs  and  particle  trajectories  can  be 
seen  further  by  comparing  the  Lagrangian  fields  to  a  tracer 
field.  Though  it  is  not  materially  conserved,  relative  humid¬ 
ity  may  be  used  as  a  tracer  field  (by  advecting  the  field  for¬ 
ward  in  time  and  pretending  there  are  no  sources  or  sinks). 
The  relative  humidity  tracer  field,  at  700  hPa  in  Fig.  8  and 
500  hPa  in  Fig.  9,  shows  the  initial  moisture  content  of  air 
parcels  entering  the  pouch  and  confirms  that  some  of  the  air 
entering  the  pouch  is  initially  dry  while  the  air  exiting  the 
pouch  is  initially  moist.  Attracting  and  repelling  LCSs  are 


overlaid  in  blue  and  red  respectively.  The  relative  humidity 
tracer  field  indicates  that  dry  air  is  already  entering  laterally, 
Fig.  8a.  From  2  September  to  4  September,  moist  air  persists 
at  the  circulation  center,  but  some  dry  air  passes  to  the  cen¬ 
ter  through  the  open  pathway.  By  3  September,  the  relative 
humidity  tracer  field,  showing  relative  humidity  values  on  1 
September,  shows  an  intrusion  of  dry  air  in  the  pouch  which 
entered  from  the  southwest  along  the  second  repelling  LCS 
branch.  The  lateral  dry  air  intrusion  is  confirmed  in  the  ana¬ 
lyzed  ozone  field  by  a  broad  intrusion  on  the  west  side  of  the 
pouch,  Fig.  8c,  d.  By  4  September,  the  air  from  the  north  has 
forced  the  attracting  LCS  southward  and  deformed  the  por¬ 
tion  of  the  repelling  LCS  wrapping  around  the  west  side  of 
the  pouch.  The  net  result  is  that  the  dry  air  intrusion  from  the 
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north  can  enter  the  pouch  directly  without  wrapping  around 
from  the  west. 

There  are  strong  similarities  between  the  tracer  field  and 
the  attracting  structures  shown  in  Fig.  8a,  b.  In  contrast,  the 
relative  humidity  analysis  data  shown  in  Fig.  7  does  not  pro¬ 
vide  much  information  about  the  transport  of  particles.  A 
similar  scenario  can  be  seen  for  pre-Karl  in  Figs.  11  and  13. 
The  intrusion  of  dry  air  continues  over  the  next  3  days,  and 
by  5  September,  the  relative  humidity  tracer  field  shows  that 
initially  dry  air  has  entered  the  pouch  while  moist  air  has  ex¬ 
ited  the  pouch  to  the  west.  The  expulsion  of  moist  air  seen 
in  the  relative  humidity  tracer  field  occurred  at  the  same  time 
that  the  repelling  LCS  emerged  at  the  pouch  center.  This  re¬ 
pelling  LCS  and  the  absence  of  an  attracting  LCS  provides  a 
mechanism  for  moist  air  to  be  expelled  from  the  center  to  the 
west  of  the  storm. 

The  RH  tracer  field  and  LCSs  can  also  be  compared  with 
the  initial  locations  of  particles  that  will  enter  the  pouch  in 
a  two  day  period.  These  particle  locations  are  computed  as 
the  final  location  of  backward  trajectories  seeded  inside  a 
3  degree  radius  circle  surrounding  the  pouch  center  from  2 
days  in  the  future.  The  backward  trajectories  for  ex-Gaston 
are  shown  in  Fig.  6  for  700  hPa  from  1  to  6  September.  Over 
this  time  period,  there  are  many  trajectories  from  the  north 
and  northeast  of  the  pouch  that  enter  the  circulation.  Com¬ 
paring  these  locations  to  the  RH  tracer  field  indicates  that 
many  of  the  trajectories  with  initially  low  relative  humidity 
values  reach  the  pouch  center. 

The  length  of  trajectory  paths  during  a  72  h  integration 
time  are  shown  in  Fig.  8e,  f.  The  gradient  of  the  curve  lengths 
aligns  with  the  attracting  LCS.  Though  the  initial  positions  of 
trajectories  on  opposite  sides  of  the  LCS  may  be  very  close, 
the  length  of  their  paths  is  very  different.  The  moist  air  from 
the  south  takes  a  very  short  path  to  the  circulation  center 
while  the  dry  air  to  the  north  takes  a  very  long  path,  with 
average  particle  speeds  of  approximately  12  m  s_1  for  a  72  h 
integration  time. 

Vertical  structure 

The  vertical  structure  of  Gaston  is  shown  on  1  and  2  Septem¬ 
ber  in  Fig.  10.  The  low-level  attracting  LCSs  (cyan)  are 
nearly  orthogonal  to  the  LCSs  at  700  hPa  and  above  (blue), 
showing  that  trajectory  paths  between  vertical  levels  can  be 
almost  orthogonal  even  when  the  Eulerian  boundaries  are 
aligned.  After  1  September,  the  attracting  LCS  is  rotated  and 
tilted,  disrupting  the  vertical  alignment  of  the  pouch  as  the 
pouch  boundary  at  700  hPa  and  above  is  south  of  the  bound¬ 
ary  below  700  hPa.  The  tilt  between  500  hPa  and  700  hPa  can 
be  seen  by  the  lack  of  alignment  of  the  eastern  hyperbolic 
trajectory,  Figs.  7c  and  8d.  A  tilted  pouch  makes  it  possible 
for  dry  air  which  has  come  from  the  north  to  enter  the  low 
level  circulation  from  aloft.  However,  the  two-dimensional 
trajectory  analysis  as  well  as  the  thermodynamic  analysis  of 
Smith  and  Montgomery  (2012)  suggest  that  the  dry  air  is 


from  the  north.  On  2  September  the  already  shallow  pouch 
deteriorates  further  as  the  repelling  LCS  at  700  hPa  near  the 
stagnation  point  detaches  from  the  LCSs  below  700  hPa,  and 
the  attracting  LCS  at  700  hPa  detaches  from  the  LCSs  below 
700  hPa.  The  attracting  LCSs  at  700  hPa  actually  align  or¬ 
thogonal  to  those  below  700  hPa.  The  repelling  LCS  on  the 
west  side  of  the  pouch  has  also  become  stronger  at  this  time 
and  intersects  the  attracting  LCS  near  the  pouch  center. 

Further  change  in  the  vertical  alignment  on  the  pouch  oc¬ 
curs  from  1  to  5  September.  While  the  pouch  boundary  re¬ 
mains  in  the  same  location  relative  to  the  pouch  center  below 
700  hPa,  the  boundary  moves  southward  at  700  hPa  as  par¬ 
ticles  from  the  north  push  the  boundary  southward.  As  the 
LCS  is  pushed  southward,  the  circulation  is  tilted,  Fig.  7. 

In  Fig.  10,  the  FTLE  ridges  are  overlaid  on  relative  hu¬ 
midity  fields  for  ex-Gaston  and  pre-Karl  at  levels  of  500  hPa, 
700  hPa,  and  925  hPa  .  The  pouch  boundary  of  ex-Gaston  is 
displaced  and  rotated  between  vertical  levels,  while  pre-Karl 
has  far  greater  vertical  coherence.  For  both  storms,  the  man¬ 
ifolds  associated  with  the  saddle  to  the  west  of  the  pouch  do 
not  coincide  with  LCSs  at  multiple  vertical  levels. 

Based  on  the  foregoing  results,  the  LCSs  are  related  to  the 
import  of  dry  air,  the  vertical  tilting  of  the  vortex,  a  constric¬ 
tion  of  the  ITCZ  moisture  source,  and  the  expulsion  of  moist 
air  when  a  repelling  LCS  is  present  at  the  center.  The  re¬ 
pelling  LCS  is  particularly  important  since  it  is  the  signature 
of  a  vortex  which  is  turned  inside-out,  shifting  from  a  pe¬ 
riod  of  net  inflow  to  net  outflow.  The  LCS  view  differs  from 
the  study  of  Schubert  et  al.  (1999),  which  showed  a  vortic- 
ity  profile  that  was  radially  inverted  due  to  shear  instabili¬ 
ties.  Because  the  radial  profile  of  vorticity  is  approximately 
monotonic,  the  profile  does  not  support  these  instabilities. 
The  thermodynamic  implications  of  a  dying  vortex  are  im¬ 
portant  also,  as  the  planar  outflow  limits  the  ability  of  con¬ 
vection  to  transport  moisture  upward.  For  a  developing  vor¬ 
tex,  the  amplification  of  vorticity  occurs  ultimately  through 
vortex  tube  stretching  associated  with  convection.  However, 
the  presence  of  a  repelling  LCS  wrapping  around  the  vortex 
core  indicates  that  convective  activity  is  diminished  since  the 
air  which  has  been  transported  inward  is  now  mixed  or  trans¬ 
ported  outward.  In  Rutherford  et  al.  (2012),  the  presence  of  a 
repelling  LCS  marks  the  end  of  the  convective  lifecycle  of  a 
three-dimensional  vortical  hot  tower.  The  splitting  of  attract¬ 
ing  LCSs  that  enclose  vortices  has  also  been  seen  in  dying 
vortices  (O’ Farrell  and  Dabiri,  2010). 

5.2  Evolution  of  Karl 

5.2.1  Summary 

Karl  was  classified  as  a  tropical  disturbance  on  8  September, 
as  an  African  Easterly  wave  merged  with  a  surface  trough 
to  the  north  of  South  America.  From  10  to  14  September, 
the  near-surface  circulation  gradually  increased  in  associa¬ 
tion  with  an  approximate  diurnal  cycle  of  convective  activity 
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Fig.  6.  The  stable  (red)  and  unstable  (blue)  manifolds  of  hyperbolic  stagnation  points  and  streamlines  (white)  of  the  co-moving  frame  are 
overlaid  on  the  relative  humidity  field  ( %)  for  Gaston  from  1  September  to  6  September  at  700  hPa.  Green  dots  indicate  particles  that  will  be 
within  3  degrees  of  the  pouch  center  in  48  h. 
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Fig.  7.  Attracting  LCSs  (blue)  and  repelling  LCSs  (red)  are  represented  as  level  contours  of  the  FTLE  field  at  80  %  of  the  maximal  FTLE 
values  and  overlaid  on  the  relative  humidity  field  ( %)  at  700  hPa  for  ex-Gaston  from  1  to  6  September. 
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Fig.  8.  Attracting  LCSs  (blue)  and  repelling  LCSs  (red)  are  represented  as  level  contours  of  the  FTLE  field  at  80  %  of  the  maximal  FTLE 
values  and  overlaid  on  the  relative  humidity  tracer  field  at  700  hPa  (a,  b)  and  the  ozone  mass  fraction  (kg  kg-1)  (c,  d)  for  ex-Gaston  on  3 
September  (a,  c)  and  5  September  (b,  d).  Backward  trajectory  lengths  are  shown  for  3  and  5  September  in  (e)  and  (f). 
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Fig.  9.  Attracting  LCSs  (blue)  and  repelling  LCSs  (red)  are  represented  as  level  contours  of  the  FTLE  field  at  80  %  of  the  maximal  FTLE 
values  and  overlaid  on  the  relative  humidity  tracer  field  at  500  hPa  (a,  b)  and  the  ozone  mass  fraction  (kg  kg-1)  (c,  d)  for  ex-Gaston  on  2 
September  (a,  c)  and  3  September  (b,  d).  Trajectory  lengths  are  shown  for  2  September  and  3  September  in  (e)  and  (f). 
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Fig.  10.  Attracting  LCSs  (blue  at  700  and  500  hPa,  cyan  at  925  hPa)  and  repelling  LCSs  (red  at  500  and  700  hPa,  magenta  at  925  hPa)  are 
represented  as  level  contours  of  the  FTLE  field  at  80  %  of  the  maximal  FTLE  values  and  overlaid  on  the  700  hPa  relative  humidity  (%)  field 
for  ex-Gaston  on  2  and  3  September  in  (a)  and  (b),  and  for  pre  Karl  on  12  and  13  September  in  (c)  and  (d). 


over  the  four  day  period  of  time  before  exhibiting  a  persistent 
region  of  central  dense  overcast  on  the  day  of  development  of 
an  organized  circulation  center,  14  September  (Montgomery 
et  al.  2011).  Karl  intensified  to  a  strong  tropical  storm  on  15 
September,  and  into  a  hurricane  on  16  September,  becoming 
a  category  3  storm  on  17  September,  shortly  before  landfall 
in  Veracruz,  Mexico.  The  progressive  increase  of  circulation 
is  shown  in  the  time  series  of  area-averaged  vorticity  within 
the  pouch  region  at  various  vertical  levels  in  Fig.  4. 

5.2.2  Lagrangian  analysis 

An  analysis  of  the  LCSs  and  their  relation  to  moisture  fields 
show  that  pre-Karl  was  far  more  favourable  for  development 
due  to  greater  vertical  alignment  of  the  LCSs,  and  a  more 
favourable  location  of  LCSs  relative  to  the  pouch  center.  The 
absence  of  a  repelling  LCS  enclosing  the  pouch  center  is  also 


important  since  there  is  no  pathway  for  particles  to  escape  the 
pouch.  The  relative  humidity  tracer  fields  show  that  the  par¬ 
ticles  entering  the  pouch  are  more  moist  than  those  entering 
Gaston. 

The  environment  of  pre-Karl  is  very  moist  at  925  hPa 
around  the  pouch  center  (not  shown),  so  that  any  entry  of 
environmental  air  into  the  pouch  will  not  reduce  the  over¬ 
all  moisture.  At  700  hPa,  Karl  is  still  very  moist,  Fig.  11, 
but  shows  some  dry  air  to  the  north  and  west  of  the  pouch. 
The  relative  humidity  structure  of  pre-Karl  at  500  hPa  on 
12  September,  Fig.  12a  is  very  similar  to  ex-Gaston  on  1 
September,  Fig.  7a. 

5.2.3  LCSs  at  700  hPa 

Figure  11a  and  b  shows  a  repelling  LCS  (red)  that  appears 
to  the  north  of  the  pouch  on  12  September,  and  is  intersected 
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by  an  attracting  LCS  on  13  September.  On  12  September, 
the  LCSs  form  an  enclosure  of  a  small  region  in  the  pouch 
center.  The  attracting  LCS  remains  to  the  north  of  the  pouch 
and  forms  a  teardrop  shaped  boundary  around  the  pouch  by 
14  September.  Evidence  of  strong  cyclonic  rotation  can  be 
seen  on  15  September  as  both  attracting  and  repelling  LCSs 
wrap  around  the  circulation  center,  Fig.  11. 

The  orientation  of  the  LCS  is  almost  90°  different  than 
for  ex-Gaston,  and  the  difference  can  most  easily  be  seen  by 
the  location  of  the  intersection  of  LCSs,  which  is  located  to 
the  northeast  of  the  pouch  for  Karl,  Fig.  lOb-d,  instead  of  to 
the  southeast  for  ex-Gaston,  Fig.  7a-d.  In  addition,  the  pouch 
boundary  of  pre-Karl  is  further  from  the  pouch  center  and  re¬ 
quires  a  longer  path  for  dry  air  to  reach  the  pouch  center.  The 
longer  pathway  can  be  seen  on  13  to  14  September,  Fig.  1  lb— 
c  as  a  dry  air  region  has  wrapped  from  the  north  around  the 
west  side  of  the  storm,  and  begins  to  enter  from  the  south¬ 
west.  The  configuration  of  the  LCSs  for  pre-Karl  still  allows 
moist  air  from  the  southeast  to  be  entrained. 

LCSs  at  500  hPa 

Figure  12  shows  the  LCS  structure  at  500  hPa  which  is  sim¬ 
ilar  to  700  hPa,  with  an  attracting  LCS  to  the  north  of  the 
pouch  center.  By  13  September,  a  repelling  LCS  intersects 
the  attracting  LCS  at  a  similar  location  as  the  intersection 
at  700  hPa.  Since  the  LCSs  align  vertically,  they  represent 
a  two-dimensional  boundary  within  the  three-dimensional 
flow  field.  The  location  of  the  attracting  LCS  to  the  north  of 
Karl  is  further  north  of  the  pouch  center  than  for  ex-Gaston, 
Fig.  12.  The  LCS  boundary  is  not  forced  southward  toward 
the  pouch  center  by  vertical  wind  shear,  so  the  boundary 
maintains  vertical  alignment. 

The  LCSs  protect  Karl  from  importing  air  from  the  north 
with  relative  humidity  levels  of  below  60  %  from  10  to  12 
September.  On  13  September,  Figs,  lib  and  12b,  the  LCSs 
at  both  500  hPa  and  700  hPa  formed  a  boundary  to  dryer  air 
wrapping  from  the  north  side  of  the  pouch  and  entering  from 
the  west.  The  attracting  LCS  of  the  eastern  hyperbolic  tra¬ 
jectory  differs  from  the  similar  LCS  seen  in  Gaston  in  that 
it  does  not  split  and  unzip.  In  addition,  the  repelling  LCS  of 
the  eastern  hyperbolic  trajectory  did  not  isolate  the  southern 
region  of  high  moisture.  Thus,  the  ITCZ  remains  intact  as  a 
potential  source  of  fuel. 

A  complete  enclosure  around  the  pouch,  surrounding  most 
of  the  3  degree  radius  circle  around  the  pouch  center  can  be 
seen  on  15  September,  Figs,  lid  and  12d,  though  a  small 
intrusion  is  later  entrained  at  both  500  and  700  hPa  through 
lobe  dynamics.  This  intrusion  can  also  be  seen  in  the  ozone 
field,  Fig.  13,  but  this  entrainment  is  much  slower  than  the 
open  pathway  that  appeared  for  ex-Gaston,  and  did  not  occur 
until  after  genesis. 


Trajectory  analysis 

The  relative  humidity  tracer  fields  at  700  hPa  and  500  hPa 
demonstrate  that  the  air  entering  the  pouch  center  of  Karl 
was  moist  at  both  levels,  Figs.  13  and  14  (left  column),  while 
dry  air  was  blocked  by  the  northern  pouch  boundary.  The 
location  of  the  boundary  to  the  north  of  the  pouch  shows  that 
moist  air  is  drawn  northward,  which  can  also  be  seen  by  the 
increase  in  the  size  of  the  region  of  low  ozone  air,  Figs.  13 
and  14  (right  column). 

The  length  of  trajectory  paths  is  very  different  for  Karl 
than  for  Gaston.  While  the  trajectories  for  Gaston  traveled 
the  longest  distances  outside  the  pouch  center,  the  trajectories 
for  Karl  that  recirculated  in  the  pouch  center  had  the  longest 
paths,  Figs.  11c,  d  and  12c,  d. 

Figure  15  shows  the  final  locations  of  backward  trajecto¬ 
ries  originating  in  a  3  degree  circle  around  the  pouch  center. 
While  many  particles  originate  outside  the  pouch  boundary, 
they  come  from  regions  of  high  relative  humidity.  In  contrast 
to  Gaston,  the  trajectories  approximately  follow  the  Eulerian 
streamlines,  though  there  are  still  particles  that  cross  Eulerian 
streamlines. 

The  pouch  averaged  Lagrangian  values  are  similar  on  all 
vertical  levels,  Fig.  5a,  c.  This  property  is  in  contrast  to  a 
disruption  at  700  hPa  and  above  for  ex-Gaston  that  causes 
lower  averaged  LCS  values.  The  average  Lagrangian  values 
increase  rapidly  as  pre-Karl  begins  to  develop  on  15  Septem¬ 
ber,  as  the  LCSs  form  the  boundaries  of  closed  circulation. 

Vertical  structure 

The  vertical  alignment  of  LCSs  is  very  different  for  pre- 
Karl  and  ex-Gaston.  Figure  10  shows  the  vertical  structure 
of  the  LCSs  for  pre-Karl  at  500,  700  and  925  hPa  in  b  and 
d.  The  attracting  (cyan)  LCSs  at  925  hPa  are  already  aligned 
with  the  attracting  LCSs  at  700  and  500  hPa  (blue)  on  12 
September  and  maintain  their  alignment  until  pre-Karl  de¬ 
velops.  In  contrast,  the  vertical  structure  of  the  LCSs  of  ex- 
Gaston,  Fig.  10a,  c,  showed  low-level  LCSs  perpendicular 
to  the  LCSs  at  700  hPa  and  above.  We  note  that  the  vertical 
alignment  of  LCSs  does  not  necessarily  require  the  vertical 
alignment  of  Eulerian  manifolds. 

6  Conclusions 

In  this  study,  we  have  introduced  finite-time  Lagrangian 
methods  for  time-dependent  flows  to  the  study  of  tropical  cy¬ 
clogenesis.  These  methods  were  applied  to  developing  tropi¬ 
cal  disturbances  located  along  the  wave  trough  axis  of  East¬ 
erly  waves.  The  cat’s  eye  boundary  of  the  circulation  was 
shown  to  be  a  robust  feature  even  in  the  time-dependent 
context  of  the  real  tropical  atmosphere.  However,  the  time- 
dependent  boundaries  of  the  cat’s  eye,  which  control  trans¬ 
port  in  and  out  of  the  cat’s  eye,  appear  very  different  in 
the  Lagrangian  fields  than  they  do  in  the  Eulerian  fields. 
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Fig.  11.  Attracting  LCSs  (blue)  and  repelling  LCSs  (red)  are  represented  as  level  contours  of  the  FTLE  field  at  80  %  of  the  maximal  FTLE 
values  and  overlaid  on  the  relative  humidity  (%)  fields  for  pre-Karl  at  700  hPa  from  12  September  to  15  September  (a-d).  The  backward 
trajectory  lengths  are  shown  on  12  September  and  13  in  (e)  and  (f). 
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Fig.  12.  Attracting  LCSs  (blue)  and  repelling  LCSs  (red)  are  represented  as  level  contours  of  the  FTLE  field  at  80  %  of  the  maximal  FTLE 
values  and  overlaid  on  the  relative  humidity  (%)  fields  for  pre-Karl  at  500  hPa  from  12  to  15  September. 


Though  the  eastern  hyperbolic  stagnation  point  and  mani¬ 
folds  of  the  co-moving  cat’s  eye  could  be  associated  with 
persistent  LCSs,  the  western  hyperbolic  stagnation  point  did 
not  provide  persistent  LCSs.  When  the  LCSs  from  the  west¬ 
ern  stagnation  point  existed,  a  boundary  was  formed  where 
transport  into  the  pouch  was  controlled  by  lobe  dynamics. 
However,  the  absence  of  LCSs  and  existence  of  only  a  single 
hyperbolic  trajectory  changed  the  Lagrangian  flow  configu¬ 
ration,  which  allowed  a  direct  pathway  for  dry  air  to  enter  the 
circulation. 

In  this  setting,  the  Lagrangian  reference  frame  has  ad¬ 
vantages  over  both  Eulerian  and  co-moving  frames.  While 
the  co-moving  frame  following  the  tropical  disturbance 
showed  closed  streamlines,  the  Eulerian  reference  frame  did 
not.  However,  streamlines  from  the  co-moving  frame  were 
crossed  by  trajectories,  as  expected.  In  contrast,  the  La¬ 
grangian  boundaries  are  based  on  particle  trajectories,  per¬ 


sistent  through  the  time-variation,  and  are  impermeable  to 
particle  trajectories. 

Attracting  and  repelling  LCSs  from  FTLE  fields  were 
found  in  both  a  non-developing  and  developing  disturbance 
to  be  the  boundaries  at  the  edge  of  the  pouch  associated  with 
the  cat’s  eye.  The  alignment  and  location  of  the  LCSs  defined 
pathways  that  determined  the  interaction  of  air  masses  with 
the  tropical  disturbance. 

6.1  Application  to  a  developing  and  non-developing 
disturbance 

The  Lagrangian  techniques  were  applied  to  a  developing  and 
non-developing  disturbance  observed  during  the  PREDICT 
experiment.  We  have  specifically  examined  the  cases  of  ex- 
Gaston,  which  did  not  redevelop,  and  pre-Karl,  which  did  de¬ 
velop  into  a  hurricane.  The  LCSs  showed  boundaries  which 
constrained  the  transport  of  moisture  and  vorticity.  There  was 
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Fig.  13.  Attracting  LCSs  (blue)  and  repelling  LCSs  (red)  are  represented  as  level  contours  of  the  FTLE  field  at  80  %  of  the  maximal  FTLE 
values  and  overlaid  on  the  RH  tracer  field  at  700  hPa  for  pre-Karl  from  13  to  15  September  (a,  c,  e).  The  LCSs  are  overlaid  on  the  ozone 
mass  fraction  (kg  kg-1)  in  (b),  (d),  and  (f). 
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Fig.  14.  Attracting  LCSs  (blue)  and  repelling  LCSs  (red)  are  represented  as  level  contours  of  the  FTLE  field  at  80  %  of  the  maximal  FTLE 
values  and  overlaid  on  the  RH  tracer  field  at  500  hPa  for  pre-Karl  from  13  September  to  15  September  (a,  c,  e).  The  LCSs  are  overlaid  on  the 
ozone  mass  fraction  (kg  kg-1)  in  (b),  (d),  and  (f). 
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Fig.  15.  The  stable  (red)  and  unstable  (blue)  manifolds  of  hyperbolic  stagnation  points  and  streamlines  (white)  of  the  co-moving  frame  are 
overlaid  on  the  relative  humidity  (%)  field  for  Karl  from  10  September  to  13  September  at  700  hPa.  Green  dots  indicate  the  locations  of 
particles  that  will  be  within  3  degrees  of  the  pouch  center  in  48  h. 


a  difference  in  LCSs  and  their  position  relative  to  moisture 
for  the  developing  and  non-developing  cases. 

Gaston  and  Karl  had  similar  structure  at  lower  levels  of 
850  hPa  and  below  but  were  far  different  at  700  hPa  and 
above.  Gaston  did  not  redevelop  due  to  the  intrusion  of  dry 
air,  the  expulsion  of  moist  air  with  high  vorticity,  and  a  re¬ 
moval  of  a  moisture  source  from  the  ITCZ.  In  addition,  the 
presence  of  a  repelling  LCS  in  the  inner  core  occurred  as 
Gaston  began  to  expel  moist  air.  The  LCSs  showed  the  pre¬ 
cise  pathways  and  boundaries  which  were  responsible  for 
all  of  these  processes.  The  vertical  coherence  of  the  pouch 
boundary  is  significantly  degraded  for  ex-Gaston.  The  pouch 
boundary  becomes  tilted  and  the  upper-level  boundary  be¬ 
comes  irreversibly  detached  due  to  impinging  vertical  shear. 

The  LCSs  for  pre-Karl  were  far  more  favourable  for  de¬ 
velopment.  The  LCSs  blocked  the  intrusion  of  dry  air,  and 
retained  a  connection  of  moist  air  inflow  from  the  ITCZ.  A 


trajectory  analysis  showed  the  role  of  LCSs  in  organizing 
transport,  and  showed  how  Karl  retained  high  moisture  at  all 
levels  while  Gaston  did  not. 

6.2  Further  use  of  Lagrangian  methods  in  cyclogenesis 
cases 

The  computation  of  manifolds  by  first  locating  the  hyper¬ 
bolic  trajectory  is  preferable  for  computing  lobe  dynamics. 
However,  the  lack  of  persistent  stagnation  points  and  am¬ 
biguity  of  proper  reference  frame  make  this  computation 
difficult.  Computing  manifolds  directly  is  difficult  to  im¬ 
plement  in  forecasting  applications.  Still,  the  approxima¬ 
tions  to  Lagrangian  manifolds  through  scalar  field  meth¬ 
ods  are  present  in  both  developing  and  non-developing  dis¬ 
turbances,  and  the  computation  of  FTLE  fields  in  an  auto¬ 
mated  manner  makes  the  use  of  these  methods  applicable  for 
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forecasting  purposes.  These  methods  are  particularly  suitable 
for  automation  since  they  require  no  knowledge  of  pouch 
translation  speed.  Though  Lagrangian  methods  are  devel¬ 
oped  for  three-dimensional  flows,  vertical  velocity  data  for 
sparse  temporal  global  model  output  would  introduce  a  high 
degree  of  uncertainty  to  trajectory  motions. 

6.3  Discussion  and  outlook 

LCSs  show  both  the  planar  pathways  for  stirring,  and  reveal 
a  complete  picture  of  the  vertical  structure  of  the  wave-pouch 
region.  While  the  planar  mixing  in  these  cases  is  supported 
by  an  analysis  of  thermodynamics  (Smith  and  Montgomery, 
2012),  further  study  of  the  three-dimensional  kinematics  and 
flow  boundaries  will  help  to  explain  the  roles  of  stirring  and 
downdrafts  for  importing  dry  air  and  expelling  moist  air.  The 
methods  used  in  this  study  will  be  applied  to  additional  pre¬ 
genesis  cases  from  the  2010  PREDICT  experiment  to  help 
identify  the  reasons  for  the  development  non-development  of 
these  pre-genesis  cases.  Real-time  analysis  using  these  meth¬ 
ods  will  be  explored  also. 
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